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Abstract 

We study the diffusive dynamics of a hard-sphere ffuid confined between parahel smooth hard 
wahs. The position-dependent diffusion coefficient normal to the walls is larger in regions of high 
local packing density. High density regions also have the largest available volume, consistent with 
the fast local diffusivity. Indeed, local and global diffusivities as a function of the Widom insertion 
probability approximately collapse onto a master curve. Parallel and average normal diffusivities 
are strongly coupled at high densities and deviate from bulk fluid behavior. 

PACS numbers: 66.10.Cg, 61.20.Ja 
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Beginning with Einstein's theory of Brownian motion, it has been reahzed that diffusion 
provides an excellent approximation to the motion of molecules in a bulk liquid for times 
much longer than intervals between molecular collisions. In molecularly confined systems, 
the situation is more complicated. On one hand, particle structures emerge that result in a 
spatially inhomogeneous density profile; and on the other hand, the relaxation time of local 
density fluctuations may become faster than the time required for particle motions to become 
"diffusive." Alternatively, long-lived correlations can prevent the system from entering into 
a diffusive regime altogether, as in single-file transport Nevertheless, in simulations of 
partially confined systems such as fluids in two-dimensional (2d) sht pores, diffusion is indeed 
observed parallel to the (quasi) infinite confining planes However, the situation is less 

clear for motions in the perpendicular direction. Even if diffusion were a useful description 
of the single-particle motions, one would expect the diffusion coefficient to be spatially 
inhomogeneous. That, combined with the spatially varying density profile and the confining 
boundaries, essentially eliminates the usual way of estimating diffusion coefficients from the 
mean square displacement as a function of time 4]- As a consequence, diffusion in highly- 
confined environments has largely remained unexplored, despite its relevance for micro- and 
nanofiuidic devices [s]. 

We use a recently proposed propagator-based formalism to estimate the position- 
dependent diffusion coefficients self-consistently from simulation trajectory data [g^. For 
diffusion, the propagator (or Green's function) G{z, At\z', 0) for single-particle displacements 
along the coordinate z normal to the confining walls is assumed to satisfy the Smoluchowski 
diffusion equation, 

f 4{Bx(.)e-«i^|e-«Gl) (1) 

with P = l/ksT, Boltzmann's constant, and T the absolute temperature. Spatial dis- 
cretization of Eq. ([T]) 0] results in a master equation that describes the single-particle dy- 
namics between neighbouring intervals along z. Local free energies F{z) and local diffusiv- 
ities D_i_{z) are then determined self-consistently from the dynamics observed in molecular 
simulations through Bayesian inference j^. 

We apply this formalism to test whether diffusion provides a quantitative description of 
the single-particle dynamics perpendicular to the confining planes for a fluid confined in a 2d 
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slit pore. We will also explore whether empirical relations between the fluid density and the 
diffusion coeflicient identifled for diffusion in bulk and parallel to the walls |3y transfer 
trivially to perpendicular diffusion. 

Hard spheres (HS) conflned between hard walls are arguably the most basic model of 
conflnement. Nevertheless, essential physics of inhomogeneous fluids is captured, such as 
pronounced local density variations [Sj or shifted fluid-solid phase boundaries with respect 
to bulk [9]. Also, the unambiguous deflnition of quantities like accessible volume, and 
rigorous ways to calculate it for hard sphere fluids [lO|, can be helpful in elucidating the 



underlying microscopic mechanisms. Moreover, theoretical and computational predictions 



for this model system can be tested in experiments on colloids 

We use discontinuous molecular dynamics (DMD) simulations to generate dynamic tra- 
jectories for our model system. To simplify the notation, dimensionless quantities will be 
reported, obtained by appropriate combinations of a characteristic length (HS particle diam- 



eter a) and time scale {a^/rri(3^ where m is the particle mass). The packing fraction = 7rp/6 
is deflned in terms of the density p based on the total (rather than center-accessible) vol- 
ume jsj. The DMD simulations each involved = 3000 identical HS particles. Periodic 
boundary conditions were applied in all directions for the bulk fluid and in the x and y 
directions for the conflned fluid. In the conflned system, perfectly reflecting, smooth hard 
walls were placed aX z = ±i7/2. The transverse self-diffusivity D\\ was obtained by fltting 
the long-time (t ^ 1) behavior of the average mean-squared displacement of the particles 
to the Einstein relation (Ar^) = 4D||t, where Ar^ corresponds to the mean-square displace- 
ment per particle in the x and y directions. To calculate D±{z) for a spatially discretized 
Eq. ([I]) 0, we use a bin size of 1/10 to divide the space in the z direction. 

To test if diffusion captures the dynamics normal to the walls, we compare the long time 
propagators G{z,t\z' of the Markovian model to the simulation data. G{z,t\z' is the 
conditional probability that a particle starting from position z' at time is found at z at 
a later time t. As shown in Fig. [T] for two different packing fractions and for a pore size 
H = 3, excellent agreement is found between the diffusive model and the MD data over six 
orders of magnitude in the propagator. 

Figure [2] compares Dj_{z) and the local density p{z) to explore the effect of wall conflne- 
ment on the local normal diffusivity. Remarkably, we flnd that the Dj_{z) is large where 
p{z) is high (except near the walls, where p{z) drops sharply). This low diffusivity near the 
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FIG. 1: The conditional probability G{z,t\z' ,0) of observing a particle at a position z at time t if 
it started at z' at t = 0. Results are shown for pore size H = 3, packing fractions (j) = 0.05 (left) 
and 0.35 (right), and times t = 1, 2, and 10 (bottom to top). The observation time t = 1 is used 
to obtain parameters of the diffusion model Eq. ([T|) (lines). Simulation results for G{z,t\z' ,0) for 
different z' are shown as symbols, where z' varies from —0.05 to —0.95 (symbol o to 0) in intervals 
of 0.1. (For reference, mean collision frequencies for the bulk hard-sphere fluid are approximately 
852 and 7369 for (p = 0.05 and 0.35, respectively.) 



walls is due to the presence of impenetrable reflective wall boundaries, limiting diffusion to 
one direction. 

At first sight, the positive correlation between high local density and faster local diffusion 
may appear counterintuitive as diffusivity for bulk fluids usually decreases with increasing 
density. To understand this unexpected behavior for inhomogeneous fluids, we turn our 
attention to the physics of layer formation in the confined environments and to why the 
fluid is structured normal to the walls in the first place. As it turns out, the confined 



fluid tends to maximize its entropy by forming these layers 



12| . Similarly, the homogeneous 
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FIG. 2: Local normal diffusivity D±{z) and density profile p{z) versus position in the z direction 
for a pore size = 5 and packing fractions = 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, and 0.40 
(top to bottom for D±{z) and bottom to top for p{z)). 

hard-sphere fluid crystallizes at high enough densities to maximize entropy [isl. [l^. We note 
that the activity C, = exp(/5/i)/A^ of an equilibrated confined system is spatially invariant, 
even though the density is not (where /z is the chemical potential and A is the thermal 
wavelength). Also, for a hard-sphere fluid, C, = p{z) / -PoCsl which means that the 
local insertion probability Po{z) (or local available volume [l5|]) is directly proportional to 
the density p{z). In other words, the counterintuitive idea that dense "layers" actually have 
more available space than the gaps between them is a consequence of this simple and exact 
relationship derived more than 40 years ago by Widom lOj. If we think of diffusion as 
particles probing their surroundings for space, then Pq should indeed be a relevant quantity. 
A "prediction" would then be that D^{z) should approximately collapse as a function of 

To test this prediction, we either need the particle insertion probability Pq as a function 
of z or the activity ^ as a function of pore size H and average pore density along with the 
density profile p{z). Using the particle insertion method can be very time consuming and 
even practically limited at high densities. Instead, we use grand canonical transition matrix 
Monte Carlo (GC-TMMC) simulations to evaluate the functional relationship between the 
activity ^ and the average pore density p for a given H. The details of the GC-TMMC 




FIG. 3: Local diffusivity D±(z) at different packing fractions as a function of the local available 
volume, measured by the local insertion probability Poiz). Filled and empty symbols correspond 
to H = 5 and 3, respectively. Data for bulk are shown by a solid line. The dashed line is a power 
law fit to the D± data. The average transverse diffusivity -Dm is also shown for reference. 



method are documented in 



16j and the specific implementation details are as in 17| 



In Fig. [3], we show D± versus Pq over a wide density range (p = 0.05 - 0.4) for two pore 
sizes {H = 3 and 5). The relationship for the bulk hard-sphere fluid and for the transverse 
diffusivity are also shown on the same plot. We find that the D± data (filled and empty 
symbols) approximately collapse to a power law form which in turn is very similar to the 
bulk relationship shown by the solid line. Note that right next to the wall, this relation 
does not hold (marked by the shaded area in Fig. [3]). However, this is expected because at 
the wall particles can only diffuse in one direction (i.e., away from the wall). Overall, the 
approximate collapse in Fig. [3] supports the idea that the local available volume, probed by 
the insertion probability, is indeed a relevant quantity for diffusion. 

We also notice in Fig. [3] that the local normal diffusivity is lower than the corresponding 
bulk value, when compared at the same value of Pq. This is expected because the presence 
of the walls directly hinders normal diffusivity in ways that are not reflected by the local 
available space. Of course, this becomes a relatively small effect at high p and is reflected 
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in Fig. [3] by a convergence of the D^, D\\, and bulk diffusivity at low Pq. The reason for the 
slightly higher Du stems from a higher chemical potential required to achieve the average 



pore density equal to the bulk and has been discussed in detail elsewhere |l7l |. 



We can also use our formalism to explore t 
and normal directions. In an earlier study 



le coupling between diffusion in the transverse 
17l |. it was observed that D\\ of the confined 



HS fluid shows pronounced negative deviations from bulk fluid behavior for relatively high 
density (e.g., = 0.4) and small pore widths {H < 3). These deviations have an oscilla- 
tory dependence on H with slower diffusion occuring for pore sizes that do not naturally 
accomodate an integer number of particle "layers" in the density profile. Based on such 
information, it was hypothesized [l7| that one might also expect a coupling between single- 
particle dynamics in directions parallel and normal to the confining walls. Here, we test 
this idea by calculating the average normal diffusivity {D±) = J^^^'^ D±pdz/ J^^'^ pdz and 
comparing it to D\\ in Fig. |H As can be seen, the if-dependent oscillation of {Dj_) closely 
tracks that of D\\. Remarkably, the oscillations in D\\ and {D±) as a function of H follow the 
pore- width dependence of the phase boundary (t>f(H) between the confined solid and fluid HS 



17]. These correlations illustrate the con- 



system joj , and also the excess entropy per particle 
nection between dynamic and thermodynamic properties arising from packing frustration. 
In particular, these results provide strong support for the picture that confinement-induced 
frustration reduces diffusive particle motions in both principal directions. 
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FIG. 4: The average normal (blue circles) and transverse (red diamonds) difFusivity versus the pore 
size H is shown for a packing fraction (p = 0.40. The corresponding bulk diffusivity is shown as 
a solid horizontal line. 
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